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BALLISTIC  RESEARCH  LABORATORIES 
REPORT  NO.  658 


Karpov/ Irr 

Aberdeen  Proving  Ground,  Md. 
18  November  1947 


THE  ACCURACY  OF  DRAG  MEASUREMENTS  AS  A  FUNCTION 
OF  NUMBER  AND  DISTRIBUTION  OF  TIMING  STATIONS 


ABSTRACT 

General  formulae  are  given  for  the  error  to  be  ex¬ 
pected  in  the  drag  coefficient  of  a  projectile  whose  flight  is 
observed  over  a  given  range  containing  n  timing  stations 
distributed  in  an  arbitrary  manner.  It  is  assumed  that  the 
time-distance  relation  can  be  represented  by  a  cubic  poly¬ 
nomial  of  the  form 

4  2  3 

t  •  a^  +  SjZ  ♦  SgZ  +  SgZ 

It  is  further  assumed  that  the  mean  errors  in  time  and  dis¬ 
tance,  e.  and  e  respectively,  are  independent  and  constant 

V  Z 

at  each  observing  station. 

Illustrative  examples  are  included. 

A  proof  of  optimum  distribution  of  timing  stations  for 
drag  determination  is  given  in  the  Appendix. 
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INTHODUCTION 

The  Spark  Range  of  the  laboratory  contains  25  spark  stations  which  are  arranged,  in  groups  of  five, 
over  a  distance  of  280  feet.  Every  odd  station  is  wired  to  record  the  time  of  passage  of  the  projectile  on 
the  drum  camera  chronograph,  13  timing  stations  in  aL.  In  addition,  stations  1,  7,  11,  13,  15,  19  and  25  are 
connected  to  six  megacycle  electronic  counters  which  record  six  time  Intervals.  The  two  timing  systems 
serve  as  a  check  on  each  other.  Distribution  of  station  is  illustrated  in  the  schematic  drawing  below.  The 
chronograph  timing  stations  are  marked  by  crosse-;  the  counter  stations  are  marked  by  dots.  In  both  timing 
systems  distribution  of  timing  stations  is  symmetric  with  respect  to  the  mid-point  of  the  range. 


With  so  many  timing  stations,  distributed  over  a  relatively  short  range,  the  drag  coefficient  can  be 
determined  with  a  high  degree  of  accuracy  with  almost  any  reasonable  symmetric  distribution  of  the  timing 
stations,  but  the  actual  distribution,  when  full  range  is  used,  appears  to  be  near  the  optimum,  i.e., leading 
to  the  mallest  error  in  the  drag. 

However,  in  actual  practice,  not  infrequently  some  timing  stations  misfire  thus  rendering  the  distri¬ 
bution  assymmetric.  Moreover,  for  special  firings  the  range  is  sometimes  shortened  by  cutting  the  last  or 
the  last  two  groups  of  stations.  Under  these  circumstances  it  is  desirable  to  know  a  priori  how  the  accuracy 
of  drag  determination  is  affected  by  using  a  shortened  range.  If  the  accuracy  is  affected  too  adversely  a  re¬ 
distribution  of  the  remaining  timing  stations  would  be  necessary. 

The  general  problem  of  finding  an  optimum  distribution  of  n  timing  statlons,for  drag  determination, 
is  of  interest  and  frequently  arises  in  practice.  Work  already  has  been  done  on  this  problem  for  cases  of 
equal  spacing  of  timing  stations,  and,  more  generally.for  symmetric  distribution  ">f  timing  stations.^  An 
analysis  of  error  in  the  drag  coefficient  for  an  arbitrary  distribution  of  the  timing  stations  is  a  simple  gen¬ 
eralization  of  the  symmetrical  case  and  is  a  straight  forward  process.  This  generalization  and  the  applica¬ 
tion  of  the  results  to  the  Spark  Range  are  the  subjects  of  this  report. 


1.  J.  L.  Kelley;  Unpublished  work.  Dr.  Kelley  considers  the  case  of  symmetric  distribution  of  timing  stations 
with  time-dfkance  relationship  given  by  a  quadratic  polynomial. 

Note:  It  is  my  understanding  that  Professor  J.  E.  McShane  has  considered  the  more  general  problem  of  the 
optimum  distribution  of  timing  stations.  Unfortunately,  his  work  has  never  been  published  and  was  not  avail¬ 
able  to  me. 
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OKRTf  ATION  OF  TBS  FORMULAS 

We  shall  assume  that  the  Ume-dlstance  relationship  can  be  represented^  by  a  polynomial.  Moreover, 

we  shall  assume  that  at  each  timing  station  the  accidental  errors  in  the  time  observations  are  characterized 

by  a  mean  error  e^,  and  the  accidental  errors  in  the  observation  of  distance  are  characterized  by  a  mean 

error  e  .  The  respective  mean  errors  are  the  same  at  each  station  and  are  independent, 
z 

We  shall  confine  our  discussion  only  to  quadratic  and  cubic  polynomials  since  there  are  both  theoreti- 

V  I 

cal  and  practical  reasons  against  using  higher  order  polynomials  for  the  spark  range.  Suppose,  therefore, 
that  the  time«distance  relation  is  given  by 

*0  ♦  V  *  • ' 

where  the  origin  of  the  z  coordinates  is  arbitrary  but  is  usually  taken  at  the  mid>point  of  the  range,  or  at  the 
average  of  the  z  coordinates. 

With  n  timing  stations  we  form  normal  equations  in  the  usual  manner,  and,  after  eliminating  a^,we 
arrive  at  the  following  set 


‘^ll*!  *  <^12*2  *  ^13*3  “  S 
^22^2  ^^23^3  •  *2 
®33*3  "  S 


•  c. 


'tj  ''Ji 


where 


2 

Cjj  ■  nj z 

C22-nX*^  -(Sz¥ 

cjj-njx* 

- 

.  XzX*^ 

Cgg-nXt®  -(X»V- 

nlzt  -  Jt  X* 

tg  -  n  X 

-  X  ^  X  tg  •  h  X  ®^t  -  X  ^  X  * 

We  shall  need  the  following  auxiliary  quantities 


^12 ’^ll  *^22*  ^^12 
^13  *  *^23  ‘  °12®13 

^23  "‘^11  S3  “  Ss* 
Cj  •^12^23’ ^13^ 


also 


^33"®13  ^22  ”‘^12  ^23 
^44  *  ^22  *^33  ■  ®23 
^2  ’^12^44*^3^ 


6 


Let  us  consider  the  timing  error  first.  After  solving  the  normal  equations  for  the  constants  a^,  a^, 
and  a^  we  precede  to  find  errors  e^.  e2.and  e^  of  these  constants  as  functions  of  e^,the  timing  error.  The 
operations  to  be  performed  are  the  following: 


where  i  •  1,  2.  3 
k-  1,  2,....n 


the  ejfpression  for  each  constant  is  partially  differentiated  with  respect  to  each  t,the  result  squared,  and  the 

2 

squared  expressions  summed  over  all  stations  •  a  standard  least  square  procedure.  The  final  results  are 
given  in  the  following  table. 

Table  I.  Error  in  constants  a^  a2,  and  a^  as  functions  of  the  timing  error  e^,  and  the  number  and  distribution 
of  timing  stations. 

A.  Arbitrary  distribution  of  stations 

2 

The  quantities  tabulated  are  (e./e,) 


Error 


Time-distance  relation 


Quadratic 


Cubic 
«®22  ^44 


"Cji  C23 


"^22  ^12 


B.  Symmetric  distribution  of  stations 

Under  this  condition  the  sums  of  all  odd  powers  of  z  vanish  and  the  above  equations  assume  much 
simpler  form. 


Error 


Quadratic 


Cubic 


Zz‘2z”.(2*T 


nJz’-fZzV 


n  2  z’  -  ( 2  z*)* 

2z^2z*-(ZaV 


2.  For  example,  see  W.  E.  Demlng:  Statistical  Adjustment  of  Data,  John  Wiley  k  Sons  Inc.,  1043. 


For  symmetric  distribution  with  equal  spacing  of  stations  the  relations  6  re  further  simplified.  Let 
the  number  of  timing  stations  be  n  •  2m  +  1,  always  an  odd  number,  the  spacing  h,  and  the  overall  range  R. 
Then  h  ■  R/2m  and  the  various  sums  of  the  coordinates  can  be  written  as: 

,4 


2z2  .  JR!  2  r^ 
(2m)'‘  r-1 


Zz 


6  2R 


(2in)' 


we  write 


C»  Symmetric  and  Equal  spacing  of  stations 


Error 


Quadratic 


2)2 


2R 

(2m? 


2  r 
1 


6  m 


2  r 

1 


6 


Cubic 


2R^  It 


(imf  2r 

2T  -  (  2>^V 


2m  -f  1 


4R^ 


the  same  as 

quadratic 


2R®  2r^2r®.(2rV 


all  summations  are  to  be  performed  from  1  to  m. 

It  Is  to  be  noted  that  error  In  the  velocity,  which  Is  proportional  to  e^,  varies  Inversely  as  the  first 
power  of  the  range;  similarly  error  In  the  drag  coefficient,'  varies  Inversely  as  the  square  of  the  range; 
and  error  In  the  derivative  of  K^,  e^,  varies  Inversely  as  the  cube  of  the  range. 


The  products 


f. 


are  Independent  of  the  range  and  depend  only  on  the  number  of  timing  stations.  These  are  tabulated  below  and 
show  the  sensitivity  of  various  errors  to  variation  In  the  number  of  timing  stations: 

Table  H.  Dep^ence  of  errors  on  the  number  of  timing  stations.  E(yial  spacing  of  stations. 


2m  4^  1 

Quadratic 

Quadratic 

and 

Cubic 

Quarttc 

and 

*1 

Cubic 

^2 

*1 

*3 

Qulntlc 

^2 

6 

1.28 

4.28 

3.80 

18.87 

7 

1.13 

3.93 

3.07 

14.70 

18.67 

0 

1.03 

3.b6 

2.71 

13.68 

14.18 

11 

0.96 

3.41 

2.48 

12.72 

12.78 

13 

0.89 

3.22 

2.37 

12.04 

11.81 
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The  error  In  Increases  nearly  fourfold  by  us\ng  quartlc  or  qulntlc  polynomials.  In  order  to  attain 
the  same  accuracy  In  drag  determination  with  quartlc  or  qulntlc  polynomials  as  with  quadratic  or  cubic  It 
would  be  necessary,  therefore,  to  double  the  range.  Increasing  the, number  of  the  timing  stations  Is  much  less 
effective  mears  of  reducing  the  error  In  ag  than  lengthening  the  range. 

Distance  error.  Distance  error  can  be  easll  ’  taken  care  of  by  the  following  device.  We  shall  assume 
tliat  distance  error  e^  can  be  expressed  as  equivalent  timing  error  by  the  relation 


where  v  Is, the  velocity.  We  already  have  C  «-  expressions  for  errors  In  the  constants  a^,  ag,  and  a^  as 
functions  of  e^^  the  true  timing  error.  The  error  In  these  constants  due  to  distance  error,  or  equivalent 
timing  error  e^g,  can,  therefore,  be  obtained  at  once  by  replacing  e^^  by  e^g.  The  total  timing  error,  therefore, 
the  two  errors  being  assumed  Independent,  Is 

J  J  J 

•t  -  •tl  ^  *12 


IMtOM  Df  ▼.  Xj,.  AMD  dX^^dv 

Neglecting  gravity,  the  equation  of  motion  of  the  projectile  can  oe  written  as 


dv  V  ^ 

Tt  ■  *  rr 


where  P.  «  ,  p being  the  air  density,  m  Is  the  mass,  and  d  Is  the  maximum  diameter  of  the  projectile. 

^  pd 

If  the  ttme>dtstance  relation  Is  given  by 

'•*o*v**/*v’ 

then  It  can  be  easily  verified  that  at  a  •  0.  the  following  relations  are  true; 


o  a. 


2a, 

*0-^1  -r 


1  ^  V. 

-PT  -JT  •=*2-  — ir- 


Let  the  percentage  mean  errors  In  velocity.  In  K^,  and  In  Us  derivative  be  Pj,  pg,  and  Pj  respective¬ 
ly.  Then  , 


Pi  -100 

Pg.lOO 


2  ^211 

te)  -(^'l 


For  all  practical  purposes,  however,  the  last  term  is  much  smaller  than  the  first  ard  with  ample  accuracy 


we  can  write 


by 


It  can  also  be  easily  shown  that  the  percentage  error  in  the  derivative  of  S'  p^,  is  adequately  ^ven 


The  expressions  for  e^,  end  e^  already  have  been  given.  These  are  functions  of  the  number  and  distri* 
button  of  the  timing  stations,  the  length  of  the  range,  and  the  total  timing  error  e^  and  can  be  computed. 

Let  a  be  the  velocity  of  sound,  M  •  ^  the  lAach  number,  then,  using  the  dafiniMons  of  constants  a^,  a2. 
and  a^  the  percentage  errors  can  be  expressed  as  follows! 

■  100  a  M  e^ 

Pg-lOO  ^  aPjasg 
6Fj  a 

h  "  ^  *3 

S  an~ 

Consider  the  percentage  error  in  K^:  it  is  proportional  to  M/K^  and  since  in  supersonic  velocities 
usually  decreases  with  increasing  M,  pg  increases  rather  rapidly  at  higher  velocities,  pg  is  also  pro¬ 
portional  to  Pj  or,  by  definition  of  F^,  it  is  proportional  to  the  mass  of  the  projectile  and  is  inversely  pro¬ 
portional  to  the  air  density,  and  to  the  square  of  the  maximum  diameter  of  the  projectile.  In  addition,  Pg  is 
proportional  to  Sg  which,  in  turn,  depends  upon  the  accuracy  of  the  instnimentation  as  represented  by  the 
timing  error  e^  and  the  distance  error  e^  and  upon  the  number  and  distribution  of  the  timing  stations. 

Finally,  Pg  depends  on  the  order  of  the  polynomial  chosen  to  represent  the  time-distance  relationship;  the 
higher  the  order  of  the  polynomial  the  greater  the  error  in  ag  (  •  Sg)  and,  hence,  in  Pg.  Therefore,  in  order 
to  maintain  the  tame  accuracy  in  with  polynomial  of  higher  order  either  the  number  of  timing  stationa 
should  bs  tncrsassd,or,  mors  sffectlvsly,  the  range  be  incressed.  As  an  tll’jstrstion,  ths  psrcentags  srror 
In  has  bean  eomputsd  for  ths  case  of  s  model  of  standard  Ibhmm  ahsll  MlOl  which  has  been  fired  in  the 
range  at  earioue  Mach  numbera.  Although  the  firinga  were  done  only  up  to  Mach  number  2.6,  the  error  curve 
was  computed  to  M  •  6.0  by  means  of  the  ereptrtcal  functtctf* 

Q  -  ^1  ♦  KpM^  -  s  ♦  bM 

For  this  mode)  s  •  0.9606  and  b  •  0.1320. 

Figure  1  shesvs  separately  the  errors  due  to  distance,  to  time,  and  the  total  error.  It  ie  apparent  that 
at  sabsottle  veloeUtes  the  distance  error  is  the  more  Important;  at  htgher  velocities  the  error  due  to  time 
predoninatea. 
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OPTIMUM  DISTRIBUnW  OF  TIMING  STATIONS 

We  have  ?een  in  the  preceding  section  that  the  error  in  is  proportional  to  the  error  in  a^,  which, 
in  turn,  with  a  given  accuracy  of  observations,  depends  only  upon  the  number  and  distribution  of  the  timing 
stations.  The  question  naturally  arises,  therefore,  whether  there  is  an  optimum  distribution  of  the  n  timing 
stations  which  will  lead  to  the  smallest  error  in  the 

For  the  case  of  symmetrical  distribution  of  the  timing  stations  when  the  time«distance  relationship 
is  given  either  by  quadratic  or  cubic  polynomials,  Dr.  H.  G.  Landau  has  shown  that  optimum  distribution  is 
attained  if  stations  are  grouped  at  each  end  of  the  range  with  one  group  in  the  middle.  In  fact,  if  the  number 
of  timing  stations  n  is  divisible  by  four,  the  best  distribution  calls  for  placing  one  half  of  the  stations  in  the 
middle  group,  with  one  quarter  of  the  stations  at  each  end.  The  proof  of  this  is  to  be  found  in  the  Appendix, 
In  general,  if  k  be  a  factor  such  that  n  •  4k  m 

where  m  can  be  either  0,  1,  2,  or  3,  and  the  range  be  two  units  lo.ng,  the  following  table  shows  the  required 
optimum  distributions. 

Table  Hi.  Optimum  distribution  of  timing  stations 
n  Number  of  stations  at 


z 

-  -  1 

z  «  0 

2  »  1 

t 

4  k 

k 

2k 

k 

4 

n 

4  k  4-  i 

k 

2k  4  1 

k 

4n 

n^  -  1 

4  k  4  2 

k 

2k  4  2 

k 

4n 

k  4  1 

2k 

k  4  1 

S  A 
n  -  4 

4  k  4  3 

k  4  1 

2k  4  1 

k  4  1 

4n 

■5 - 

n  -  1 

To  quote  from  the  Appendix:  "It  should  be  pointed  out  that  it  will  not  be  possible  to  fit  a  cubic  in  z 
with  exactly  these  spaclngs  of  stations,  since  only  three  values  of  z  are  given.  Because  it  is  physically  im¬ 
possible  to  put  more  than  one  station  at  one  position,  the  above  spacing  can  only  be  approximated.  .  .  .". 

In  the  Spark  Range  the  stations  are  arranged  in  groups  of  five,  so  we  can  use  11  timing  stations  for 
which  the  theoretical  optimum  distribution  requires  placing  3  stations  at  each  end  of  the  range  and  5  in  the 
middle.  The  following  table  compares  values  for  the  theoretical  optimum, for  physically  achievable 
distribution  in  the  range  with  stations  five  feet  apart,for  our  usual  distribution  of  every  odd  station  being  a 
timing  station,  and  also  using  seven  electronic  counter  stations; 


No.  of  Timing 
Stations 

4  2 

10^  X  — i 
«t 

Optimum 

11 

.295 

Ftosslble  in  the  Range 

11 

.305 

Usual  distribution  (drum  ca,mera) 

13 

.364 

Electronic  counters 

7 

.461 
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U  the  above  figures  are  translated  Into  the  percentage  error  In  K^,  Pg,  of  model  155mm  shell  MlOl  for 
example,  the  effect  on  Pg  of  various  distribution  appears  as  follows: 

Model  155mm  shell  MlOl  Pg  at  M  2.0 


Optimum 

0.109 

Possible  In  the  Range 

0.113 

Usual  distribution  (drum  camera) 

0.135 

Electronic  counters 

0.171 

Thus  theoretical  error,  using  only  the  electronic  counters  3tatlons,ls  50%  larger  than  the  optimum  error 
attainable  In  the  range.  However,  for  various  practical  reasons  U  Is  undesirable  to  segregate  the  timing 
stations  as  Indicated  by  the  optimum  distribution.  Moreover,  the  whole  error  Is  so  small  that  even  seven 
counters  provide  ample  accuracy. 

The  following  additional  few  simple  examples  Illustrate  the  effect  of  various  distributions  on  Sg.  Sup¬ 
pose  the  length  of  the  range  Is  10  units.  The  table  below  gives  various  distributions  and  corresponding  values 
of  Sg/Sj. 


Distribution 

No.  of  stations 

V*t 

•6,  “2  0,  0,  0,  +2^,  +6 

7 

.0960 

-6,  C,  0,  0,  +5 

5 

.0366 

•6,  -1,  0,  +1,  +6 

6 

.0376 

•6,  -2,  0,  •*■2,  +6 

5 

.0405 

-6,  -3,  0,  +3,  ■♦■5 

5 

.0453 

-5,  0,  0.  +6 

4 

.0400 

One  station  misfires 

-5,  X,  0,  +2,  +5 

4 

.0434  Quadratic 
,0490  Cubic 

The  table  shows  that  U  Is  possible  to  attalp  as  good  an  accuracy  with  fewer  stations  properly  located 
.  as  with  greater  number  placed  less  judiciously. 


THE  ORDER  OF  TIME -DISTANCE  POLYNOMIAL 


The  order  of  the  polynomial  representing  the  time-distance  relationship  in  a  given  range  should  be 
such  that  the  next  higher  term  should  be  less  than  the  error  in  time  measurement.  This  is  a  necessary 
condition.  Tnus  If 

.  2  3  m 

t  -  a  +  0,2  +  a„z  ♦  a„z  +  .  .  .  +  a  z 
0  1  z  d  m 

•6 

and  the  time  is  measured  tc  the  accuracy  of  10  seconds,  the  above  condition  requires  that 


m 


a  z 
m 


<  10 


-6 


ntax 
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If  the  z’s  are  reckon««l  from  the  mid  point  of  the  range,  z  ”  r  .  therefore 

m&x  Ct 


Given  vs  M,  the  various  coefficients  can  be  computed.  Thus 

/  'S 


M 


24F“  a  L  dM‘ 


Ir  (• 


In  the  following  table  a^  and  a^  are  tabulated  for  various  Mach  numbers  for  the  model  155mm  shell  MIDI 
whose  vs  M  graph  Is  given  in  Figure  2. 


M 

^3 

"4 

1.5 

4.582x10"^^ 

-3.60x10 

2.5 

1.688 

-0.156 

3.5 

0.790 

+0.032 

4.5 

0.439 

+0.034 

The  amd  its  derivatives  were  computed  from  the  Q  function.  The  Fj  for  this  shell  is  750  feet.  For  esti¬ 
mating  the  order  of  the  polynomial  the  above  table  can  be  used  with  sufficient  accuracy  for  other  shell 
provided  the  tabular  values  are  multiplied  by  the  ratio  of  750  to  F^  of  the  new  shell,  raised  to  appropriate 
power  as  indicated  in  the  formulae.  ^  ^ 

In  the  following  table  and  a^  are  tabulated  for  the  spark  range,  =*  280  feet,  and 

the  transonic  range,  Rg  -  700  feet.  The  transonic  range  of  the  laboratory,  which  is  being  built,  will  contain 
25  photographic  stations,  arranged  in  groups  of  fives,  over  a  distance  of  700  feet.  For  this  example  it  was 
assumed  that  3-inch  shells  were  to  be  fired  in  this  range  for  which  Fj  »  2660  feet. 


M 

^4(r 

R  «  280 
155mm  model 

R  -  700 

3" 

R  -  280 

R  -700 

1,5 

12,55x10'® 

15.59x10'® 

-13.80x10"® 

-4.84x10 

2.5 

4.63 

5.75 

-  0.60 

-0,21 

3.5 

2.16 

2.68 

+  0.12 

+0.04 

4.5 

1.20 

1.49 

+  0.13 

+0.05 

On  the  basis  of  aforementioned  criteria,  therefore,  in  both  spark  and  transonic  ranges,  the  fourth  power  term 
can  be  safely  omitted.  ^ 


It  is  to  be  noted  that  although  the  cubic  term  is  retained,  the  determination  of 


nevertheless,  is 


very  poor.  This  can  be  seen  from  the  values  of  Pg,  the  percentage  error,  tabulated  on  page  15  for  the  model 
155mm  shell. 
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P3 

1.5 

71.2 

2.5 

133.8 

3.5 

261.5 

4.5 

457.9 

Clearly,  the  derivative  cannot  be  accurately  determined  from  a  single  round. 

It  is  of  interest  to  see  how  well  the  time-distance  relationship  could  be  approximated  by  a  polynomial 
if  there  were  no  observational  errors.  Advantage  can  be  taken  of  the  observed  fact  that  in  supersonic  range 
of  velocities,  variation  of  the  drag  coefficient  with  ^4ach  number  is  accurately  represented  by  the  Q  function. 

3 

Neglecting  gravity,  therefore,  the  equation  of  motion  can  be  integrated  explicitly.  The  result  is: 

-Ct  2 

z  -  -  At  +  Bln  ( ^^2 - iJ  \ 

V  D  -  1  ' 

where  A,  B,  C,  and  D  are  functions  of  M  at  zaO  and  of  the  slope  and  intercept  of  the  Q  function  of  a  particular 
type  of  projectile. 

The  computations  were  done  for  the  cone-cylinder  model,  round  1730,  which  was  fired  at  The 

z  coordinates  were  computed  by  the  above  fcrmula  from  the  observed  times  thus  freeing  the  z,  t  values  from 
observational  errors.  Polynomials  of  various  degrees  were  fitted  to  these  values  by  least  squares  with 
origin  of  the  coordinates  being  kept  always  at  station  13,  the  midpoint  of  the  range;  for  both  symmetric  and 
asymmetric  alstribution  of  timing  stations.  The  following  table  contains,  for  each  polynomial,  the  resulting 
K^,  Its  percentage  error  P2  computed  from  the  residuals  in  the  usual  manner,  and  the  mean  error  e^,  in 
microseconds,  of  an  equation  of  unit  weight.  The  quality  of  fit  can  be  judged  from  the  size  of  e^. 

Table  IV 

Quality  of  fit  of  z,  t  values  by  a  polynomial 


Distrib. 

No.  of 
stations 

Quadratic 

Cubic 

Quartic 

Symmetric 

7 

.1250 

.1250 

.1227 

P2 

.25 

.14 

.81 

®o 

1.77 

l.Ol 

1.05 

Asymmetric 

6 

.1240 

.1240 

.1249 

stations  1-19 

P2 

« 

0 

.14 

.95 

®o 

.20 

.24 

.53 

Asymmetric 
stations  7-25 

0 

.1256 

..235 

.1165 

P2 

.40 

.30 

5.9 

e 

0 

1.48 

.66 

3.06 

3.  I  am  indebted  to  Dr.  A.  C.  Charters  for  this  formula. 
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On  the  whole,  at  this  Mach  number,  the  cubic  seem  to  give  the  best  .  however,  quadratic  could 
have  been  used  equally  well.  Experience  shows,  however,  that  over  the  whole  range  of  Mach  numbers, 
and  especially  through  the  transonic  range,  cubic  representation  leads  to  more  consistent  results.  It  is  to  be 
noted  that  using  the  observed  values  of  (z,t)  the  resuduals  in  time,  and  hence  the  mean  error,  will  usually  be 
augmented  by  Inlierent  inadequacy  of  the  polynomial  to  fit  the  (z,  t)  function  as  shown  in  Table  IV.  Thus 
with  the  observed  (z,  t)  values,  using  symmetric  distribution  and  cubic  polynomial,  the  results  are: 

Kj^  .1251 
P2  .29 

e  2.06 

o 

with  the  mean  error,  e^,  twice  as  large.  However,  one  microsecond  of  this  error  is  due  to  failure  of  the 
polynomial  to  fit  the  data,  and  only  the  remainder,  therefore,  is  attributable  to  errors  in  time  and  distance. 


USE  OF  POLYNOMIAL  REPRESENTATION  IN  THE  TRANSONIC  REGION 

Few  remarks  should  be  made  concerning  the  use  of  polynomial  fitting  of  (z,  t)  data  in  the  transonic 
region  where  varies  rapidly  with  the  Mach  number. 

It  should  be  recalled  that  a  cubic  representation  of  the  (z,  t)  data  takes  care  of  the  linear  part  of  the 
variation  of  by  the  cubic  term.  In  the  spark  range,  the  retardation  of  projectiles  over  280  feet  seldom 
exceeds  .08  Mach  numbers  and  is  usually  considerably  less,  so  the  linear  approximation  of  variation  of 
vs  M  is  adequately  taken  care  of  by  the  cubic  term. 

Another  alternative  sometimes  is  being  advocated,  namely,  first  to  differentiate  once,  numerically,  the 
(z,  t)  data,  and  represent  the  velocities  by  a  quadratic  in  z.  From  the  equation  of  motion,  again  neglecting 


gravity,  we  have 


at  z  -  C 


The  alx)ve  procedure,  could  perhaps  be  used  successfully  for  projectiles  with  large  retardation  such 
as  spheres  or  irregular  fragments.  For  ordinary  projectiles  the  reduction  of  data  by  the  above  method  leads 
to  somewhat  larger  errors  in  than  str  light  polynomial  reduction. 

The  two  methixis  designated  I  and  11  respectively  have  been  applied  to  model  flrlr.gs  of  155mm  shell 
MIDI  with  results  tabulated  on  page  17.  The  values  of  with  corresponding  percentage  errors  computed 
from  residuals  are  tabulated  for  various  M. 
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Table  V 

Representation  of  z,  t  data  by  polynomial  (I),  and  by  method  II. 
lS5mm  shell  MlOl 


Rd. 

M 

Method  I 

Method  II 

Kd 

P2 

^2 

1211 

.797 

.0690 

.26 

.0682 

.30 

1394 

.934 

.0656 

.43 

.0684 

4.16 

1207 

.936 

.0709 

.10 

.0709 

.84 

1208 

.966 

.0962 

.18 

.0926 

2.01 

1389 

.993 

.1299 

.29 

.1281 

1.44 

1206 

1.017 

.1444 

.11 

.1444 

.91 

1206 

1.066 

.1698 

.18 

.1681 

1.43 

1381 

1.066 

.1640 

.19 

.1637 

1.01 

1382 

1.071 

.1646 

.09 

.1636 

.68 

1318 

1.966 

.1187 

.19 

.1171 

1.96 

1312 

1.999 

.1164 

.28 

.1206 

2.42 

Although  the  differences  in  K^’s  computed  by  two  methods  are  not  exessive  and  are  nonsystematie 
in  character, the  percertage  errors  by  the  second  method  appear  to  be  about  10  times  larger.  Representatlot 
of  (z,  t)  data  by  a  cubic  polynomial  even  in  the  transonic  range  must  be  considered  quite  Mtisfactory. 

I  wish  to  acknowledge  my  indebtedness  to  Mr.  E.  Dearden  and  Mr.  K.  Q.  Tadman  of  the  British  Branch 
for  Theoretical  Research,  Fort  Halstead,  Kent,  England. 


B.  G.  Karpov 
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APPENDIX 

OPTHSUM  DnTRlBUnON  op  tdiinq  stations 
FOR  DRAG  DETERMINATIONS 


This  appendix  by  H.  O.  Landau  gives  the  solution  of  the  following  problem;  To  determine  the  drag 

coefficient  for  a  projectile  at  the  center  of  the  Aerodynamics  Range,  n  stations  are  placed  symmetrically 

about  ihe  banter  of  the  range  and  the  position,  z,  of  each  stations,  and  time,  t,  when  the  projectile  passes 

is  detsrmined.  A  quadratic  or  cubic  in  z  is  fitted  to  t  by  least  squares,  and  the  error  in  drag  coefficient 

2 

will  be  a  minimum  when  the  error  in  the  coefficient  of  z  is  a  minimum.  How  should  the  positions  of  the 
stations  be'chosen  so  as  to  minimize  this  error? 

2 

Since  the  square  of  the  error  in  a^,  the  coefficient  of  z  ,is  proportional  to  the  reciprocal  of 


the  problem  is  to  maximize 


n 


n  Z 
i-1 


2 


8  •  n 


P 


where  the  i|  are  .the  distances  of  the  stations  from  the  center  of  the  rangs  and  the  seals  is  chosen  so  that 
the  ends  of  the  range  are  a  •  1  and  a  «  -I. 


0  <  Xj «  1  ,  (2) 


and  from  the  symmetry,  every  value  of  x^  ^  0  must  occur  an  even  number  of  times. 

8  is.  of  course,  the  second  moment  of  the  x^  about  their  mean  ttimes  the  constant,  n).  The  problem 
is  to  find  the  distribution  of  the  Xj  which  maximises  8  subject  to  (2).  ntysically,  we  must  find  the  positions 
of  n  point-masses  on  a  weightless  bar  of  unit  length  which  give  the  largest  moment  of  inertia  about  the  center 
of  gravity.  The  answer  iS  almost  obvious  intuitively:  half  the  x^  must  have  the  value  0  and  the  other  half  the 
value  1. 

The  proof  can  be  given  on  the  basis  of  the  identity. 


S  •  n 


(3) 
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which  can  be  seen  to  be  true  by  expanding 


n  A  i  ”  a  ' 

.  .  Z  X,*  .  I  .  ) 

M  '  M 

Let  the  x^  be  numbered  la  order  of  Inereastnc  macnltude 


‘x.  *  1. 


aod  euppoee  first  that  n  la  even.  Then  we  can  see  that  any  distribution  of  the  x^  fives  8  a  amaller  value  than 
that  flven  by  the  distribution  with 

Xj  ■  0  for  *  ^  j  » 


X, -Ifor  t  *  5  ♦  1 

For  any  I,  the  sum  of  the  terms  in  S  which  include  x^  are  from  (3) 


S,  -  X  (*,  -  *,  ) 

‘  j.l  ‘  ' 


We  show  that  if  i  ^  then  is  Increased  by  putttnf  x^  •>  0.  Using  (4). 

2  "  2 
8,  J  0  -  1)  x  *  ♦  X  (t  *  xT 
*  '  * 


n  »  n  j 

X  -  X  (X.  -  x^)  ♦  x^ 


"  a  “  r  2' 

-  X  (»,  •  *,)  ♦  2x,  X  (x,  -  x^)  ♦  (n  -  t  -  1)  X, 

J-t>l  *  ’  J-l*l  ^ 

V 


the  ssoowd  term  here  te  positive  and  the  last  term  la  gre^r  than  (I  •  11  x^*  because  i  ^  j 

n  «  n  n  n 

Z  X,  »  Z  (X,-X)*.0.1).^‘  . 

H.I  '  )-l.l  ' 


So  that 
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which,  using  (6)  gives,  ,  , 

n  n  ^  0 

s^^  2  X  ^  Z  X  +  2  X 

the  value  of  for  ■  0, 

In  exactly  the  same  way  It  can  be  seen  that  If  t  ^  ^  f  1,  is  increased  by  putting  >  1. 

Hence  it  follows  that  the  maximum  value  of  S  is  given  by  the  distribution  (5). 

If  n  is  a  multiple  of  4,  n  >  4k,  this  gives  the  answer  immediately:  half  of  the  have  the  value  0, 
one-qti^er  will  be  -fl  and  one  quarter  will  be  -1. 

The  modification  for  n  not  a  multiple  of  four  can  be  seen  as  follows: 

a.  n  oddc  Just  as  before  we  must  have  x^  -  0  for  I  ^  ,  and  x^  -  1  for  I  ?  -y-'  .  The  value 

of  x^j  is  not  yet  determined. 

2  The  sum  of  the  terms  which  Include  x  ,  is 


S  ,  -  T  x^  ♦  (1  -  X 

n+1  2  n+l  n+1 


Vl’  ^ 

L  T  T  J 

The  equality  is  reached  for  either  x^^^ «  0  or  1. 

■f 

For  n  ■  4k  ♦  1,  the  symmetry  condition  requires  putting  x^^^  »  0,  and  for  n  •  4k  3 
symmetry  requires  x^^j  •  1. 


b.  n  •  4k  ♦  I.  In  this  case  the  distribution  (6)  would  not  satisfy  the  symmetry  condition.  We  know 


that  for  S  to  be  a  maximum  we  must  have 


X,  -  0  for  i  ^  ,  and  x.  -  1  for  i  ^ 


and  symmetry  requires  x  -  x  „  ,  Then  just  as  in  (a)  abovt ,  it  follows  that  these  two  x’s  either  both  have 

n  n+4 

5  T 

^e  value  0  or  both  equal  1 . 

Tbe  results  are  summarizec  in  tl  e  following  table. 


Number  of  values  of  z  at 
n  z«-l  z-C  z«l  S 


4k 

k 

2k 

k 

n 

4 

4k  >  1 

k 

2k  4  1 

k 

2  , 
n  -  1 

4k  >  2 

k 

2k  4  2 

k 

2  4 

n  -  1 

k  ♦  1 

2k 

k  4  1 

4k  ♦  3 

k  ♦  1 

2k  ♦  1 

k  4  I 

n^  -  1 

4 
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It  should  be  pointed  out  that  It  will  not  be  possible  to  fit  a  cubic  In  z  with  exactly  these  spaclngs  of 
stations,  since  only  three  values  of  z  are  ^ven.  Because  It  Is  physically  impossible  to  put  more  than  one 

3 

station  at  one  position,  the  above  spacing  can  only  be  approximated,  and  the  coefficient  of  z  can  then  be  de 
termlned,  but  Us  error  will  be  large.  If  this  coefficient  Is  desired  with  any  accuracy  the  problem  must  be 
reconsidered,  taking  this  Into  account. 
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